Statistical Properties of E(5) and X(5) Symmetries 


Jing Shu!, Hong-bo Jia!, Yu-xin Liul?34* 


! Department of Physics, Peking University, Beijing 100871, China 
? The Key Laboratory of Heavy Ion Physics of the Chinese Ministry of Education, 
Peking University, Beijing 100871, China 
? Institute of Theoretical Physics, Academia Sinica, Beijing 100080, China 
^ Center of Theoretical Nuclear Physics, National Laboratory of 


Heavy Ion Accelerator, Lanzhou 730000, China 


April 14, 2017 


Abstract 


We study the energy level statistics of the states in E(5) and X(5) dynamical 
symmetries. The calculated results indicate that the statistics of E(5) symmetry 
is regular and follows Poisson statistics, while that of X(5) symmetry involves two 
maxima in the nearest neighbor level spacing distribution P(s) and the A statistics 
follows the GOE statistics. It provides an evidence that the X(5) symmetry is at 


the critical point exhibiting competing degrees of freedom. 
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The nature of “shape” phase transitions in finite quantal system is a fundamental 
issue and has been the subject of many investigations. In the last years, there have been 
many preeminent works concerning this subject in nucleus, including both theoretical 
calculations[l] P] B, H] and experimental results, f]. Particularly, two symmetries at 
the critical points of shape phase transition, called E(5) and X(5), have recently been 
proposed[] B]. It is a major breakthrough in the study of critical point behavior of 
nucleus undergoing a shape phase transition because it gives us simple analytic solutions 
while historically we have to resort to numerical calculations to describe such nucleus. 

Empirical examples of such kind of symmetries have been found soon after they 
were put forward. For instances, the !?Ba, !?pq, !5Pqd, Ru isotopes are proposed 
to be the ones in E(5) symmetry B, [10, and the Sm, Ng, 194Mo in X(5) 
symmetry [13] [[5]. These examples show that the obtained new symmetries represent 
a helpful theoretical tool which could well describe the structure of the realistic nucleus 
at the critical point of a shape transition, in a sense complementing the role played by 
the dynamical symmetries of the Interacting Boson Model(IBM)[18, [I]. However, there 
remain some interesting aspects on the statistical properties of the E(5) and X(5) sym- 
metries. (1) Whether the nuclear system displays ordered or chaotic spectrum in the new 
symmetries? (2) Does the statistics give a strong evidence to support that the obtained 
new symmetries can well describe the nucleus at the critical point in a shape transition? 
(3) The new symmetries are based on the particular potential amenable to analytic de- 
scriptions which well approximate to the "true" potential. Does the simplification of the 
potential eliminates part of the statistical properties of the realistic nucleus at the critical 
point? Or more concretely, are the results consistent with the statistical properties ob- 
tained by Alhassid and collaborates|[3, [9| in the framework of IBM? Therefore it is 


crucially important to obtain the statistics of the new symmetries and answer the above 


opening questions. In this sense, we explore the statistical properties of E(5) and X(5) 
symmetries in the paper. 


We start from the Bohr Hamiltonian[1] 


ge dp eee, "xum ev 3». Q) 
2B-0^08 OB T Riess s "Ay sin^(y — £n) id 


In the E(5) symmetry, the potential V(6,7) depends only on the deformation parameter 
B and it can be written as V(3,y) = U(B). Considering the case, in which the potential 


is simplified as an infinite well, one has the eigenvalues[/] 


Eg, = i az 


where B is a parameter, xz, is the £th zero of the Bessel function J,(z), and the order 


Ds aa (2) 


of the Bessel function is v = T + 3 with 7 being the quantum number similar to the 
irreducible representation of O(5) group in the IBM, w is the width of the well. 

In the X(5) symmetry, the potential was supposed to be a square well in the variable 
B and a harmonic oscillator in y with no £ — y couplings, which is a good approximation 


in the U(5)-SU(3) shape transition. In this case, the eigenvalues are given as|§ 


E(s, L, n4, K, M) = Eg + E}, (3) 

where 

Xs,L 2 
E a 4 
p= (2) (4) 
and 
3a 4 (K/2)? 

E= (n, 4- 1) (0/2) ' (5) 


(B2) 3 (8?) 
The x, is the sth zero of Bessel function J,(z) of (irrational) order v = [J(J + 1)/3 + 
9/4]'/2, L is the total angular momentum (with projections K on the symmetry axis and 


M on the quantization axis in the laboratory frame). a is a parameter related to the 


interaction strength of the potential U, = (3a)?y?/2 and n, is the number of 4-vibration 
quanta, while Sw is the same as that in E(5) symmetry and the (8?) is the average 
of 8? over the wave functions. Combining all variables, one obtains the most general 


expression [|B] 


E(s, L,n4, K, M) = Eo + B(x: L) + An, + CR’, (6) 
where 
A= _ 3a (7) 
(8?) 
and 


B- CE (8) 


Before analyzing the energy level statistics, we estimate the relative value of the param- 
eters A/B in the light of the “true” potential given in the Fig.1. of Ref.[B]. We can see 
from the figure that the width of the infinite well to approximate the “true” potential is 
Bw ~ 0.8, and most of the quantum states can thus be restricted in Ug < 0.5 (arbitrary 
unit). If we assume that the total energy of the quantum states are conserved in E — 0.5, 
the quantum states of the y degree are then mainly restricted in U, < 0.5, as the potential 
is completely decoupled into two separate parts in 6 and y respectively. Since the X(5) 
symmetry is obtained around y = Of], the corresponding restriction in y by the total 
energy must be very small. Taking this small value of y as 0.1, we then work out that the 
strength of the two dimensional oscillator (3a) must be more than 10 according to the 
formulation U., = (3a)?7?/2. With Eqs.(7) and (8), we obtain that the parameter A must 
be more than 25 if we take V (2) as 0.4 and the parameter B — 1.56. Consequently we 
get the relative value A/B > 16. 

Our estimate is consistent with the recent empirical example !''^Mo[L5] of X(5) sym- 
metry, where the nucleus presents the X(5) pattern not only in its ground state band, but 


also in its low-lying n. = 1,2; K = 2n, bands. From the level scheme given in the Fig.2. 
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of Ref.[I5], we can obtain the empirical values of the parameters, A ~ 800 and B ~ 30 
(in keV). The relative value is clearly in the region of A/B > 16 as we have proposed. 
Such a consistence supports that the assumption to obtain X(5) symmetry around y = 0 
in the Ref.[§] is quite reasonable since in fact, most of the quantum states in the realistic 
nucleus at the critical point of the transition from spherical to axially-deformed shape are 
restricted in y < 0.1. 

To analyze the energy level statistics in the dynamical symmetries, we take the fol- 
lowing process. At first with Eqs.(2) and (6), we calculate the energy levels of the nucleus 
in E(5) and X(5) symmetries, respectively. The parameter € and 7 in E(5) symmetry 
are truncated manually at a maximal number €( and 7”), We also truncated s and 
n, at s™ and nim) respectively in X(5) symmetry. After the spectrum (£;) has been 
determined, it is necessary to separate the smoothed average part whose behavior is 
nonuniversal and can not be described by random-matrix theory(RMT) [RJ]. To do so we 
count the number of the levels below E so that one can define a staircase function N(E) 


of the spectrum (see for example Ref. f18]) 
N(E) = Nal E) + Nga (E). (10) 


Then we fix the N,,(E;) semiclassically by taking a smooth polynomial function to fit the 


staircase function N(E). We obtain finally the unfolded spectrum with the mapping 
{Ei} = NS). (11) 


This unfolded level sequence LE;) is obviously dimensionless and has a constant average 
spacing of 1, but the actual spacings exhibit frequently strong fluctuation. 

We have used two statistical measures to determine the fluctuation properties of 
the unfolded levels: the nearest neighbor level spacings distribution(NSD) P(S) and the 


spectral rigidity A3(L). The nearest neighbor level spacing is defined as S; = (E;,1) — (E;). 
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The distribution P(S) is defined as that P(S)dS is the probability for the S; to lie within 
the infinitesimal interval [S, S + dS]. For a regular system, it is expected to behave like 


the Poisson statistics 

P(S)— e$, (12) 
whereas if the system is chaotic, one expects to obtain the Wigner distribution 

P(S) = (1/2)S exp(—1S?/4), (13) 


which is consistent with the GOE statistics[D2, D4] D3]. With the Brody parameter w in 


the Brody distribution 


P,(S) = a(1 + w)S"exp(—a 0*9), (14) 
where 
o =T[(2+w)/(1+w)]/? (15) 


and [|z] is the [ function, the transition from regularity to chaos can be measured with 
the Brody parameter w quantitatively. It is evident that w = 1 corresponds to the GOE 
distribution, while w = 0 to the Poisson distribution. A value 0 < w « 1 means an 
interplay between the regular and the chaotic. 
As to the spectral rigidity A3(L), it is defined as 
1 pL 

A3(L) = (minant [ING — Az — Bus) : (16) 
where N(x) is the staircase function of a unfolded spectrum in the interval [- 7/2, 2]. 
The minimum is taken with respect to the parameters A and B. The average denoted 
by (---) is taken over a suitable energy interval over x. Thus from this definition A3(L) 


is the local average least square deviation of the staircase function N(x) from the best 
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fitting straight line. For the GOE the expected value of A3(Z) can only be evaluated 


numerically, but it approaches the value 

A3(L) & (InL — 0.0687) /1? (17) 
for large L. And for Poisson statistics 

Ag(L) = L/15. (18) 


In our calculation of E(5) symmetry, all the degenerate states caused by additional 
quantum number are considered just as one single state[D5]. That is because we manu- 
ally introduced the additional quantum number to distinguish the degenerate states with 
the same quantum number 7. The degenerate states with different additional quantum 
numbers are statistically uncorrelated since no interaction is involved to link such states. 
If we do not remove the degeneracy, we would take the mixed ensembles (with different 
additional quantum numbers which are good quantum numbers in the dynamical symme- 
tries) into consideration and obtain the over-Poisson distribution in practical calculations 
since nearly 1/4 of the spacings of all are zero. In order to evaluate the fluctuations of the 
energy levels in one pure ensemble, these large amount of the statistically uncorrelated 
states should be removed. 

At first, we analyze the energy levels given in Eqs.(2) for the E(5) dynamical symme- 
try. The numerical results of the NSD and A; statistics show that the statistical feature 
does not depend on the angular momentum obviously. We then illustrate only those of 
the states with low spin-parity J^ = 0+ but different sets of truncated numbers € and 
7 of E(5) symmetry in Fig.1 (a), (b) and (c). Meanwhile the Brody parameter w of 
the NSD[4] is also evaluated and listed in the figures. The numerical results show that 
the statistics of the E(5) symmetry is quite close to Poisson statistics. When the number 


of the energy levels involved in the statistics (denoted by €™ and 77?) increases, the 
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statistics becomes almost exactly the Poisson-type. It indicates that the E(5) symmetry 
is a dynamical symmetry which is completely integrable in classical limit. 

We then analyze the statistics in the X(5) symmetry. We can easily find that the 
energy eigenvalues in the X(5) symmetry contain two terms, one is the square of the zeros 
of the Bessel function in 8, another is the solution of two dimensional oscillator function 
in y. In order to display how the relative strength of the two terms affects the energy 
level statistics, we calculate the statistics in different values of the relative strength A/B, 
the results are displayed in Fig.2. All of the relative strength A/B satisfy the condition 
A/B > 16 as we estimated above. It is apparent that when the parameter A/B is not very 
large, two maxima[R5] appear in the NSD statistics (see the left part of Fig.2. (a), (b)), 
which is rather different from the fluctuation properties in other cases except for those in 
the U(5) dynamical symmetry with collective backbending[25, BG]. In theoretical point 
of view, the X(5) symmetry describes the critical point nuclei in a spherical to axially 
deformed shape phase transition analytically. The structure of its energy levels must be 
not the same everywhere. In fact, the term B(X. L)” is approximately a rotational term 
of s, whereas the term An, is a vibrational term of n,. The simultaneous appearance of 
these two terms induces a competition and makes the ensemble with the same angular 
momentum L and its projection on the symmetry axis K involves two sequences, one of 
which is in vibration, another is in rotation. When we unfold the above spectrum {F£;(L)}, 
the two sequences are normalized with an unique total average spacing, and their maxima 
in the spacing distributions are no longer the same, consequently two maxima emerge 
because of the uniform unfolding procedure. The same mechanism and results can be 
found in U(5) dynamical symmetry with collective backbending, where the the structure 
of the yrast band changes from the U(5) vibrational states to the rotational states with 


full d-boson number ng = N, which implies that a transition of collective motion mode 


may happen as the angular momentum increases. In this statistical point of view, the 
X(5) symmetry describes well the critical point behavior of nucleus undergoing a shape 
transition from U(5) dynamical symmetry to SU(3) dynamical symmetry which involves 
competing degrees of freedom. 

In the A; statistics when A is not very large (see the right part of Fig.2. (a), (b)), it 
is clearly that the statistics follows the GOE statistics, which suggests that the competing 
strengthes in the 6 and y degrees can cause the onset of chaos in X(5) symmetry. The 
case that the degree of chaoticity reaches its maximal when different competing strengthes 
coexist and are comparable can be also found in other works (see for example Ref. B3}, R7] ). 

Looking over Fig.2.(c) and (d), one can recognize that when we increase the strength 
of the oscillator in y degree of freedom, which corresponds to the parameter A, the 
two maxima in the NSD statistics soon vanishes and both the NSD and As statistics 
become close to the Poisson statistics. It indicates that only the vibrational mode plays 
a dominant role in the system and the competition mentioned above does not exist. If we 
go to extremes to the case when the vibrational strength 3a — oo, the y degree of freedom 
no longer exists. The X(5) symmetry becomes the E(5) symmetry in a four-dimensional 
space, and its statistics is certainly regular. 

The results obtained above are consistent with Alhassid and collaborates’ work[19]. 
The E(5) symmetry lies in the transitional region between a vibration to a y-unstable 
rotation, which is regular both classically and quantum mechanically in the work. It has 
been known that the Hamiltonian in the consistent Q formalism[8] of the IBM can be 
given as 


H = eia — &Q^(x) - Q'(x), (9) 


where 


Q? (x) = (std + dts)? + x(did)?. 


By changing c/& under restriction x = — 7/2, one can describe the U(5)-SU(3) transition 
with this Hamiltonian. To discuss the shape phase transition and relate it with X(5) 
symmetry, we reparametrize the ration e/& as (1 — Q)/C. From the Ref.[l] we know that 
the flat potential near the critical point could be obtained when the parameter ¢ €(0.025, 
0.029) [I] in case that the boson number N = 10. We then transfer the parameter ¢ into 
parameter 7 used in Alhassid and collaborates’ work[I9]. The relation is 


L6 


ToT a 


and the value of 7 near the critical point is 0.770 ~ 0.796. It is clearly that the critical 
point which corresponds to the X(5) symmetry lies in the most chaotic region(see the 
Fig.1. in Ref[19]). 

Finally, we turn to the question we have raised in the beginning. In X(5) symmetry the 
nucleus exhibits chaotic spectrum, which is probably caused by the comparable competing 
strengthes in the £ and y degrees. The two maxima appearing in the NSD statistics 
suggest that the spectrum splits into two different sequences, one is rotational and the 
other is vibrational. Such a behavior illustrates that the nucleus in X(5) symmetry is at the 
critical point of the transition from spherical to axially-deformed shape with competing 
degrees of freedom. While the statistics in E(5) symmetry provides no evidence for such 
a drastic competition in the critical shape. That might be because the E(5) symmetry is 
a critical symmetry corresponding to an isolated point of second order phase transition. 
Very recently, it has been pointed out that such a second-order phase transition located 
at the triple-point of the three phases[9]. The shape coexistence region shrinks[I] as the 
parameter x in the Casten’s triangle changes from —/7/2 to 0. In deed, as a general rule, 
phase coexistence can be found only in first-order phase transitions[B0]. There is no shape 
coexistence in E(5) symmetry and the spectrum does not split into two sequences, then 


no competition exists and its statistics is uniform(no two maxima in the NSD statistics). 
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The common subalgebra O(5) in the U(5)—O(6) transition makes the statistics regular 
in this region and the ordered spectrum of nucleus in E(5) symmetry confirms the past 
work[I9]. It is worth to mention that the bifurcation effect(number of maxima more than 
one) in the NSD statistics might be a useful tool to determine whether the system is 
at the critical point of a first order shape transition or shape coexistence which involves 
competing degrees of freedom. 

In summary, we have calculated the statistics of E(5) and X(5) symmetries. The 
statistics of X(5) symmetry involves two maxima in the NSD statistics and the A statis- 
tics follows the GOE statistics, which provides an evidence that the X(5) symmetry is at 
the critical point exhibits competing degrees of freedom. The statistics of E(5) symmetry 
is regular. That is probably caused by the common subgroup O(5) in the U(5)—O(6) tran- 
sition. Since the E(5) symmetry corresponds to the second order shape phase transition 


with no shape coexistence, then no such competition exists as in the X(5) symmetry. 
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Figure 1: Comparison of energy level statistics of the states J^ = 0* in E(5) symmetry 
with different numbers of energy levels: (a) for 7? = 20, €°” = 20, (b) for 7°” = 40, 
€(™ = 40, and (c) for 7°” = 80, £? = 80. In all figures, the solid lines and dashed lines 


describe the GOE and Poisson statistics, respectively. 
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Figure 2: Comparison of energy level statistics of the states J7 = 8* and K = 6 in X(5) 
symmetry with the same manually truncated quantum number s" = 30 and nim) = 30, 
but different relative strength A/B: (a) for A/B = 20, (b) for A/B = 30, (c) for A/B = 
60, and (d) for A/B = 100. 
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